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1 We consider the time evolution of massless gravitino perturbations in Schwarzschild black 

' holes, and show that as in the case of fields of other values of spin, the evolution comes in 

£-"■■ three stages, after an initial outburst as a first stage, we observe the damped oscillations 

■ characteristic of the quasinormal ringing stage, followed by long time tails. Using the 

' sixth order WKB method and Prony fitting of time domain data we determine the 

' quasinormal frequencies. There is a good correspondence between the results obtained 

by the above two methods, and we obtain a considerable improvement with respect to 
the previously obtained third order WKB results. We also show that the response of a 
black hole depends crucially on the spin class of the perturbing field: the quality factor 
becomes a decreasing funetion of the spin for boson perturbations , whereas the oppositc 
Situation appears for fermion ones. 
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1. Introduction 



The study of the evolution of small perturbations in black hole backgrounds is a very 
interesting subject. Actually, we know that this evolution, at intermediate times, is 
dominated by damped Single frequency oscillations. These characteristic oscillations 
have been termed quasinormal modes and the associated frequencies quasinormal 
frequencies, and depend only on the Paramete rs cha racterizing the black hole, as it 
mass, electric charge and angular momentum l l | 2 | 3 [ In this sense we can say that 
black holes have a characteristic sound, resembling for example the familiär sound 
produced bay the ringing of a bell or the strum of a guitar. 

Once one realizes the importance of black holes in fundamental phy sics, o ne can 
grasp the meaning and significance of their characteristic vibrations Black 
holes have been called the hydrogen atom of general relativity, perfect comparison 
because like the hydrogen atom in quantum mechanics,a black hole, as a Solution 
of the Einstein's field equations has all the general relativistic properties embodied 
in it, but still is simple enough to be a model for starting a complete understanding 
of all the physics that go with Einstein's gravitation theory. 



1 
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There are different contexts in which the study of quasinormal modes of black 
holes app ears to be motivated: the estimation of astrophysical black hole param- 
eters 12141 j the study of the stability of these Solutions under small perturbations, 
semiclassical ways to quantize the black hole area and the estimation of thermaliza- 
tion timescales in connection with the ADS/CFT correspondence, in which a large 
static black hole in asymptotically AdS spacetime corresponds to a thermal State 
in a Conformal Field Theory ( CFT ) at its boundary, and the decay of the test 
held in the black hole spacetime corresponds to the decay of the perturbed State in 
the CFT. In this sense, we can determine the timescale for the return to thermal 
equilibrium studying its dynamics in AdS spactime, and then translating it onto 
the CFT, using the AdS/CFT conjecture^El 

Since the discovery of quasinormal oscillations of black holes, many studies have 
been done on quasinormal modes of various spin fields and a considerable vari- 
ety of analytica l, semi-analytical and numerical m ethods have been developed to 
dctcrminc t hcm '7 |SI3ll()lllll2ll3 ' 14ll{ > llfill7lJ ™. 

An interesting problem to consider is the study of quasinormal modes of s pin 
s = 3/2 fields, called Rarita-Schwinger fields. In a previous paper Shu and Shen 1 ^ 
determined the quasinormal modes of this fields in Schwarzschild spacetime using 
the third order WKB method developed by Iyer and Will 

In this paper we investigate for the first time the complete time evolution of 
gravitino perturbations in the Schwarzschild background and determine with better 
accuracy the quasinormal frequencies by Utting time domain data by a superposition 
of damping exponents, using the Prony method. We also extend previous results 
by Shu and Zheng bey ond th e third order WKB, using the sixth order approach 
developed by Konoplya The results obtained by this two approaches match 

very well, ensuring the validity of the determined quasinormal frequencies. We also 
studied the complete time evolution of scalar, electromagnetic, gravitational and 
Dirac perturbations in the mentioned spacetime, and compute numerically some 
quasinormal frequencies by the two methods mentioned above. 

The structure of this paper is as follows. In Section II we use the Newmann- 
Penrose formalism to obtain the equations for gravitino test perturbations in 
Schwarzschild spacetime, and use Separation of variables to obtain decoupled equa- 
tions for the radial and angular variables involved. The radial equation is written in 
the familiär way useful to study the quasinormal spectrum of the perturbations, ob- 
taining an analytical expression for the effective potential related to it. Section III is 
devoted to determine the time domain evolution and to compute some quasinormal 
frequencies of gravitino perturbations, and to compare them with similar calcula- 
tions for other boson and fermion fields. The last section contains our conclusions 
and presents some lines for future research on this subject. 
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2. Rarita-Schwinger equation in Schwarzschild spacetime 

In supergravity theory the gravitational Held is coupled to the massless spin 3/2 
Rarita-Schwinger field, that acts as a source of torsion and curvature. if we denote by 
^ abc ^ ne spin 3/2 field, then the supergravity Heids equations are invariant undcr 
supergravity transformations 4>abö ~~ ^ ^ 'abc + ^ BC aA > wnere is an arbitrary 
spinor field. 

When the Rarita-Schwinger field vanishes, the supergravity equations becomes 
Einstein's vacuum equations. Therefore, since the torsion and the curvature pro- 
duced by the Rarita-Schwingcr field dopend quadratically on this field, then in the 
linear approximation about a Solution with i^abc = 0' ^ ne supergravity field equa- 
tions reduces to the Ra rita- Sc hwinger field equations together with the vacuum 
Einstein's field equations I 22 | 23 | 

The Rarita-Schwinger equations for the field iPabc m a curve d background can 
be written as 

HABC = H(ABC), Habc = ® (-0 
where we have defined H A nn = V, h £ V' 4 ^a a and H A ■■ = V D r -ib A tr ^-, (the 

\ts ^ C)D BC (ß \D\C) y 

parenthesis denotes symmetrization on the indices enclosed and the indices be- 
tween bars are excluded from the symmetrization). It is possible to show that in 
an algebraically special vacuum spacetime, the contraction of H AB( j defined above, 
with a multiple principal spinor of the conformal curvature satisfies a decoupled 
equation ^H. 

Using the Newmann-Penrose notation, we can show that in a frame such that 
^2 is the only nonvanishing component of the Weyl spinor, the Rarita-Schwinger 
reduces to only two equations for the componcnts Hqqq = fii and Hm = f7 3 , 
that can be solved by Separation of variables in all the type D vacuum backgrounds. 
In terms of the differential Operators A = (D — 2e + e* — 3p — p*) (A — 37 + /Lt), 
B = (ö - 2ß - a* - 3r + tt*) (5* - 3a + tt), C = (Ä+2j-'y*+3ß+ß*) (D + 3e - p) 
and D = (S* + 2a + ß* + 3n — r*) (8 + 3ß — r) that contains the usual directional 
derivatives of the Newmann-Penrose formalism, we can put the Rarita-Schwingcr 
equations in the above mentioned frame as 

(l-fl-* 2 )fi| =0 (2) 



c-D-y 2 )n_3 =0 (3) 

We are interested in the specific case of Schwarzschild background, where the metric 
can be expressed as 

ds * = ^dt 2 - ^dr 2 - r 2 (dO 2 + sin 2 9d0 2 ) (4) 

r A 

with A = r 2 — rrn and rjj = 2M defines the location of the event horizon. The 
separable Solution of @ and ([3]) are given by 

^1 = i?i(r)S*3(6i)e i( " t+m<#,) (5) 
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Q_3 = ^R_3{r)S_s(6)e l(ult+m ' t,) (6) 

2 2v2r 3 2 2 

whcrc uj is thc frequency, m is a half integer, and the funetions R(r) and S(9) satisfy 
the ordinary differential equations 

(aS_^®1 - 4iwr) Aii?3 (r) = \A?R%(r) 

(A^l^o + AiLür) R_s.(r) = \R_s(r) (7) 

V 2 / 2 2 

£li£-§£§(0) = -A5|(fl) 
£_i£L35_3(ö) = -A5_3(ö) (8) 

2 2 2 2 

where A is a Separation constant that takes the values (£+3)(£ + l) with £ = 2, 3, 4.... 
In the above equations we have used the Operators defined by 

2)„ =d r + — + n— (InA), = d r — — h n— (In A) (9) 

£„ = de + m esc ö + n cot 6>, £) n = de — m esc ö + n cot 6> (10) 

Now introducing the tortoise coordinate r* defined by dr* — (r 2 / A)dr and making 
use of the variable Y(r) = r~ 2 R_z(r) we can put the second equation in ([7]) in the 
Standard form 

A 2 Y + PA+Y - QY = (11) 
where A± = -£-^ ± iuj, A 2 = A + A_ = A_ A + = -£j + u 2 and 

P(r) = (4A + rr H - 2r 2 ) (12) 
2r 4 

Q(r) = ^((^ + 3)(^ + l)r + r H ) (13) 

As explained in Reference it suffices to consider the equation for R_s only, and 
an entirely similar reduetion can be made with the equation for i?| . Finally. It is 
possible to reduce (fTTj) to a, Schrodinger-type equation in thc form 

^Z(r)+[u J 2 -V(r)}Z(r)=0 (14) 

where the effective potential for the Rarita-Schwinger field is given by 

V(r) = Q(r)-^ (15) 



with the funetion T(r) defined by T(r) = -pj^(£ + 3)(£ + 1)^(^ + 3)^+1) + 1 + 

£QnF(r)), and F(r) = ^Q(r). 

In Figure (|TJ) we show the effective potential for gravitino perturbations with 
two different multipole numbers t, for a Schwarzschild background with ru = 2. As 
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Fig. 1. Effective potential of gravitino modes for Schwarzschild black hole with t = 3 (top) and 
l = 4(bottom) for r H = 2. 

we can see, V has the form of a definite positive potential barrier, i.e, it is a well 
behaved function that goes to zero at spatial infinity and gets a maximum value 
near the event horizon. 

3. Time domain evolution of perturbations and quasinormal 
frequencies 

The most direct approach to study the gravitino perturbations is to solve numeri- 
cally the evolution equation associated with (fl"4"|). that is 



where the function §(t,r) results from the factorization of f2_| in ^ as fl_3 = 
&(t, r)S_ 3 (9)e tm ^. In order to integrate numerically the equation ([Tüll we use the 

technique developed by Gundlach, Price and Pulling and the result can be 
observed as the time-domain profile showed in Figure ([2]). 

As we can see, the temporal evolution of gravitino perturbation can be divided 
in three stages. The first depends on the initial conditions and the point where 
we observe the profile. After the initial outburst at the first stage, we observe the 
cxponential damping of the perturbations called quasinormal ringing, that can be 
split to the superposition of exponentially damping oscillations, represented by a 
set of complex frequencies, called quasinormal frequencies, whose real parts describe 
the actual frequency of the oscillation, while the imaginary part is the damping rate. 




(16) 
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The quasinormal modes are Solutions of the wave equation (|14[) with the specific 
boundary conditions requiring pure out-going waves at spatial infinity and pure in- 
coming waves on the event horizon. Thus no waves come from infinity or the event 
horizon. As we can see from the figure, the quasinormal ringing stage is followed by 
the so-called tails at asymptotically late times. Thus, the time evolution of gravitino 
perturbations in Schwarzschild spacetime follows the same stages as well known 
boson and fermion perturbations. To show this, we also integrated numerically the 



5 




10 20 50 100 200 



Fig. 2. Logaritmic plot of the time-domain evolution of 1 = 2 gravitino perturbations of gravitino 
modes for Schwarzschild black hole at r = llr^f. In the figure the time is measured in units of 
the horizon radious: r = t/rn- 

perturbation equations corresponding to scalar, electromagnetic, gravitational and 
Dirac perturbations in this spacetime. The resulting profiles are observed in Figure 
©. Note that in all of the cases considered, the time evolution is similar, then the 
gravitino is not the exception. In order to evaluate the quasinormal modes we used 
two different methods. The first is a semianalytical method to solve equation (TT4")) 
with the required boundary conditions, based in a WKB-type approximation, that 
can give accurate values of the lowest ( that is longer lived ) quasinormal frequencies, 
and was used in sev eral papers for the determination of quasinormal frequencies in 
a variety of Systems ^3. 

The WKB technique was applied up to first order to finding quasinormal modes 
for the first time by Shutz and Will ^"1 Latter this approach was extended to the 
third order beyon d the e ikonal approximation by Iyer and Will and to the sixth 
order by Konoplya ^^ l We use in our numerical calculation of quasinormal modes 
this sixth order WKB expansion. The sixth order WKB expansion gives a relative 
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Fig. 3. Logaritmic plot of the time-domain evolution of £ = 2 scalar (left top), electromagnetic 
(right top), dirac (left bottom) and gravitational (right bottom) perturbations in Schwarzschild 
black hole at r = llrfj. In the figure the time is measured in units of the horizon radious: r = t/rn- 



error which is about two Order less than the third WKB order, and allows us to 
determine the quasinormal frequencies through the formula 

2' 



i^—^l-^^ =n+ l -, n = 0,1,2,... (17) 
V ~ 2V o 



3=1 

where Vq, is the value of the potential at its maximum as a funetion of the tortoise 
coordinate, and V represents the second derivative of the potential with respect to 
the tortoise coordinate at its peak. The correction terms Ilj depend on the value of 
the eff ective potential and its derivatives ( up to the 2i-th order) in the maximum, 
see Eül and re ferences therein. 

The second metho d that we used to find the quasinormal frequencies was the 
Prony method HSHül f or Atting the time domain profile data by superposition of 
damping exponents in the form 



p 



k=l 



C k e~ lWkt 



(18) 



Assuming that the quasinormal ringing stage begins at t = and ends at t — Nh, 
where N > 2p — 1, then the expession (jT8]) is satisfied for each value in the time 
profile data 



ip (nh) 



p 

E 

fc=i 



C k e 



-iwkfih 



(19) 
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From the above expression, we can determine, as we know h, the quasinormal fre- 
quencies oji once we have determined Zi as functions of x n . The Prony method allows 
to find the Zi as roots of the polinomial function A(z) defined as 
v p 

A (z) = I] ( z - z k) = E a ^ zP ' m > «o = l (20) 

k— 1 m— 

It is possible to show that the unknown coefficients a m of the polinomial function 
A(z) satisfy 

p 

Solving the N — p+1 > p linear equations (1211) for a m we can determine numerically 
the roots z a and then the quasinormal frequencies. 

It is important to mention the fact that with the Prony method we can obtain 
very accurate results for the quasinormal frequencies, but the practical application 
of the method is limited because we need to know with precision the duration of the 
quasinormal ringing epoch. As this stage is not a precisely defined time interval, in 
practice, it is difficult to determine when the quasinormal ringing begins. Therefore, 
we are able to calculate with high accuracy only two or, sometimes three dominant 
frequencies. 

In Table [T] we show the results for the evaluation of the first four fundamental 
quasinormal modes for gravitino perturbations in Schwarzschild black holes, using 
the two methods mentioned above. 

The parameter entering in the calculation is the black hole mass M , that we take 
to be a unit mass. We also show for comparison the results obtained using a third 
order WKB method, that were obtained in reference With sixth order WKB 
approximation we obtain an improve of the results for the quasinormal frequencies, 
as we can easily see after comparison with the numerical results obtaining using the 
Prony fitting of the time domain data. 

The obtained results are showed in Figure For a given angular number, 
higher overtones becomes less oscilatory, i.e, the real part of the quasinormal fre- 
quencies decreases whereas the imaginary part increases. In this Situation, the qual- 
ity factor of modes, defined as the ratio T = |i?e(w)|/|7m(u;)|, decreases. On the 
contrary, modes with higher multipole number £ and the same overtone number have 
higher frequencies and the damping shows only little increments, with an effective 
increase of the quality factors. 

To test the accuracy of the Prony method to finding quasinormal frequencies by 
fitting time domain data, we also calculated the quasinormal frequencies of scalar, 
electromagnetic, gravitational and dirac perturbations in the Schwarzschild space- 
time. 

To comparative purposes, we also calculated the above freque ncie s using sixth 
order WKB method, and correct the results presented in reference for the quasi- 
normal frequencies of Dirac and gravitino perturbations, where apparently the au- 
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Table 1. Gravitino quasinormal frequencies in Schwarzschild black hole, with M = 1 
( third and sixth order WKB approximation and Prony fitting of time domain data). 



e 


n 


Third order WKB 


Sixth order WKB 


Prony 


2 





0.7346 - 0.0949« 


0.7347 - 0.0949« 


0.7348 - 0.0949« 


2 


1 


0.7206 - 0.2870« 


0,7210 - 0.2869« 


0, 7210 - 0.2869« 


3 





0.9343 - 0.0954« 


0.9344 - 0.0954« 


0.9344 - 0.0953« 


3 


1 


0.9233 - 0.2876« 


0.9235 - 0.2876« 


0.9237 - 0.2875« 


3 


2 


0.9031 - 0.4835« 


0.9026 - 0.4840« 




4 





1.1315 - 0.0956« 


1.1315 - 0.0956« 


1.1317 - 0.0956« 


4 


1 


1.224 - 0.2879« 


1.1225 - 0.2879« 


1.1227 - 0.2878« 


4 


2 


1.1053 - 0.4828« 


1.1050 - 0.4831« 




4 


3 


1.0817-0.6812« 


1.0798 - 0.6830« 




5 





1.3273 - 0.0958« 


1.3273 - 0.0958« 


1.3276 - 0.0957« 


5 


1 


1.3196 - 0.2881« 


1.3196 - 0.2881« 


1.3199 - 0.2879« 


5 


2 


1.3048 - 0.4824« 


1.3045 - 0.4826« 




5 


3 


1.2839 - 0.6795« 


1.2826 - 0.6805« 




5 


4 


1.2582 - 0.8794« 


1.2548 - 0.8832« 






Fig. 4. Gravitino quasinormal modes. 



thor uses an improper generalization of thc formula for the effective potential due 
to boson fields in Schwarzshild spacetime to the case of fermion fields. The equation 
for the radial part of boson perturbations have the form (fT4"| with the effective 
potential: 



Vi(r) = 4 



r 



(22) 
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where A = £(i+l) and the quantity ß is related to the spin-weight s of the perturbing 
field by ß — 1 — s 2 . For scalar, electromagnetic and axial gravitational perturbations 
ß has the integer value 1, 0, -3, respectively. 



Table 2. Scalar quasinormal frequencies in Schwarzschild black hole, with M = 1 
( numerical and sixth order WKB approximation and Prony fitting of tirae doraain 
data. 



1 


n 


Leaver results 


Sixth order WKB 


Prony 


2 





0.4836 - 0.0968« 


0.4836 - 0.0968« 


0.4836 - 0.0968« 


2 


1 


0.4639 - 0.2956« 


0,4638 - 0.2956« 


0.4639 - 0.2956« 


3 







0.6754 - 0.0965« 


0.6754 - 0.0965« 


3 


1 




0.6607 - 0.2923i 


0.6609 - 0.2922« 


3 


2 




0.6336 - 0.4960« 




4 







0.8674 - 0.0964i 


0.8674 - 0.0964« 


4 


1 




0.8558 - 0.2909i 


0.8558 - 0.2909« 


4 


2 




0.8337 - 0.4903« 




4 


3 




0.8032 - 0.6975i 




5 







1.0596 - 0.0963i 


1.0596 - 0.0963« 


5 


1 




1.0500 - 0.2901« 


1.0500 - 0.2901« 


5 


2 




1.0315 - 0.4873i 




5 


3 




1.0052 - 0.6899« 




5 


4 




0.9728 - 0.8995i 





The polar perturbations corresponding to a gravitational perturbing field are 
governed by the so called Zerilli potential, that is different from (|2"2"|) , but as was 
rigorously shown by Chandrasekar in the quasinormal frequencies belonging to 
axial and polar perturbations are identical. 

For the case of Dirac perturbing field, we use the form of the effective potential 
obtained by Cho in reference ^Sl, given by: 

V h M = f£ + W 2 (23) 

with 

W(r)=K^ (24) 

where the quantity k = £ + 1. Tables [5] to [5] shows the results obtained. We also 
show the third order WKB results for all the perturbations. In the case of scalar, 
electromagnetic and gravitational perturbations, we listed the numerical results 
given in the paper of Leaver 

The third order WKB results showed for the case of Dirac perturbations are a 
refinement of the Cho results ^Sl, to include one more decimal place (see also^H). 

As we can see, the Prony method gives values that are in perfect agreement 
with previous numerical results existing for scalar, electromagnetic and gravitational 
perturbations. Also we obtain a improvement of the third order WKB results using 
the sixth order semianalytic formula. 
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Fig. 5. Dependance of the quality factor upon multipole and overtone numbers for gravitino (left) 
and gravitational (right) perturbations. 



Fig. 6. Dependance of the quality factor upon spin for boson(left) and fermion(right) perturba- 
tions. 



Table 3. Electromagnetic quasinormal frequencies in Schwarzschild black hole, 
with M = 1 ( numerical and sixth order WKB approximation and Prony fitting of 
Urne domain data). 



1 


n 


Leaver results 


Sixth order WKB 


Prony 


2 





0.4576 - 0.0950« 


0.4576 - 0.0950« 


0.4576 - 0.0950« 


2 


1 


0.4365 - 0.2907« 


0,4365 - 0.2907« 


0.4365 - 0.2907« 


3 





0.6569 - 0.0956« 


0.6569 - 0.0956« 


0.6569 - 0.0956« 


3 


1 


0.6417-0.2897« 


0.6417 - 0.2897« 


0.6417-0.2897« 


3 


2 


0.6138 - 0.4921« 


0.6138 - 0.4921« 




4 







0.8531 - 0.0959« 


0.8531 - 0.0959« 


4 


1 




0.8413 - 0.2893« 


0.8414 - 0.2893« 


4 


2 




0.8187 - 0.4878« 




4 


3 




0.7876 - 0.6942« 




5 







1.0479 - 0.0960« 


1.0479 - 0.0959« 


5 


1 




1.0382 - 0.2891« 


1.0384 - 0.2890« 


5 


2 




1.0194 - 0.4856« 




5 


3 




0.9928 - 0.6876« 




5 


4 




0.9510 - 0.8967« 
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Table 4. Dirac quasinormal frequencies in Schwarzschild black hole, with M = 1 (nu- 
merical and sixth order WKB approximation and Prony fitting of Urne domain data). 



1 


n 


Third order WKB 


Sixth order WKB 


Prony 


2 





0.5737 - 0.0963« 


0.5741 - 0.0963« 


0.5743 - 0.0962« 


2 


1 


0.5562 - 0.2930« 


0.5570 - 0.2927« 


0.5571 - 0.2927« 


3 





0.7672 - 0.0963« 


0.7673 - 0.0963« 


0.7675 - 0.0962«' 


3 


1 


0.7540 - 0.2910« 


0.7543 - 0.2910« 


0.7543 - 0.2910« 


3 


2 


0.7304 - 0.4909« 


0.7297 - 0.4919« 




4 





0.9602 - 0.0963« 


0.9603 - 0.0962« 


0.9605 - 0.0961« 


4 


1 


0.9496 - 0.2902« 


0.9498 - 0.2901« 


0.9499 - 0.2901« 


4 


2 


0.9300 - 0.4876« 


0.9295 - 0.4881« 




4 


3 


0.9036 - 0.6892« 


0.9011 - 0.6925« 




5 







1.1531 - 0.0962« 


1.1531 - 0.0962« 


5 


1 




1.1442 - 0.2897« 


1.1443 - 0.2897« 


5 


2 




1.1271 - 0.4860« 




5 


3 




1.1025 - 0.6869« 




5 


4 




1.0718 - 0.8939« 





Table 5. Gravitational quasinormal frequencies in Schwarzschild black hole, with 
M = 1 ( numerical and sixth order WKB approximation and Prony fitting of time 
domain data). 



i 


n 


Leaver results 


Sixth order WKB 


Prony 


2 





0.3737 - 0.0890« 


0.3736 - 0.0889« 


0.3737 - 0.0890« 


2 


1 


0.3467 - 0.2739« 


0,3463 - 0.2735« 


0.3467 - 0.2739« 


3 





0.5994 - 0.0927« 


0.5994 - 0.0927« 


0.5994 - 0.0927« 


3 


1 


0.5826 - 0.2813« 


0.5826 - 0.2813« 


0.5826 - 0.2813« 


3 


2 


0.5517-0.4791« 


0.5516 - 0.4790« 




4 





0.8092 - 0.0942« 


0.8092 - 0.0942« 


0.8092 - 0.0942« 


4 


1 


0.7966 - 0.2843« 


0.7966 - 0.2843« 


0.7966 - 0.2843« 


4 


2 


0.7727-0.4799« 


0.7727 - 0.4799« 




4 


3 


0.7398 - 0.6839« 


0.7397 - 0.6839« 




5 







1.0123 - 0.0949« 


1.0123 - 0.0949« 


5 


1 




1.0022 - 0.2858« 


1.0022 - 0.2858« 


5 


2 




0.9827 - 0.4803« 




5 


3 




0.9550 - 0.6805« 




5 


4 




0.9208 - 0.8881« 





The dependance of the quality factor with the multipole and overtone number 
is similar as for the case of gravitino test field perturbations. In Figure ([5]) we show 
this dependance for gravitino and gravitational perturbations. For a fixed multi- 
pole number, higher overtones shows a decreasing in the quality factor, and the 
modes are more damping. In contrast, increasing the multipole number for a given 
overtone, the quality factor increases due to the increasing of the real oscillating 
frequency in a more pronounced way that the damping rate. At this point it is 
interesting to note that fermion perturbations in the black hole background shows 
higher quality factors that boson perturbations, as Figure (j6]) shows. As the spin 
weight of a boson perturbation increases, the quality factor decreases, and the black 
holes perturbed by this Heids becomes more poor oscillators. Then, the lowest qual- 
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ity factors corresponds to gravitational perturbations, and the higher to scalar test 
perturbations. 

The Situation for fermion fields is opposite. In this case the quality factor in- 
creases with the increment of the field spin, and as a consequence the highest quality 
factor belongs to the gravitino perturbations. Then, as a general fact, black holes 
perturbed by fermion fields are better oscillators than those perturbed by boson 
fields. 

4. Concluding remarks 

We have studied the evolution of Rarita- Schwinger field perturbations in a 
Schwarzschild background. Solving numerically the time evolution equation for this 
perturbations, we find similar time domain profiles as in the case of fields of other 
spins: the usual three stages in the time evolution dominated at intermediary times 
by quasinormal ringing. We determined the quasinormal frequencies by two different 
approaches, 6th order WKB and time domain Integration with Prony Atting of the 
numerical data. The two methods give close values of QNMs for well pronounced 
potential barriers. 

Also we apply the Prony method to numerically calculated the quasinormal fre- 
quencies of scalar, electromagnetic, gravitational and Dirac perturbations in the 
Schwarzschild background, with results in perfect agreement with previous numer- 
ical calculations by other authors. 

The time domain evolution of gravitino perturbations are similar to that corre- 
sponding to other fields of different spin. The behaviour of the quasinormal frequen- 
cies with respect to the black holes masses is also similar for all spin perturbations. 
However, an interesting difference occur in the oscillating behaviour of black holes 
perturbed by fields of different spin classes: if the perturbing field has integer spin, 
then the System becomes a poor oscillator as the spin increases. The opposite Situa- 
tion occurs in the presence of fermion perturbations, black holes perturbed by higher 
spin fields are better oscillators. In general, the highest quality factor is associated 
with gravitino perturbations and the smallest with gravitational ones. 

There are extensions of this work that are interesting to consider, first the gen- 
eralization to the case of charged black holes, and second, the determination of the 
changes of the gravitino quasinormal spectrum for semiclassical Solutions, due to 
the vaeuum polarization of quantized fields. The Solution of the above problems will 
be presented in future reports. 
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